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Abstract. A comparison is made between bispectral operator pairs and dual pairs of 
isomonodromic deformation equations. Through examples, it is shown how operators 
belonging to rank one bispectral algebras may be viewed equivalently as defining 1- 
parameter families of rational first order differential operators with matricial coefficients 
on the Riemann sphere, whose monodromy is trivial. By interchanging the roles of the 
two variables entering in the bispectral pair, a second 1-parameter family of operators 
with trivial monodromy is obtained, which may be viewed as the dual isomonodromic 
deformation system. 



1. Bispectral Operators. 

1.1 Bispectral Pairs. 

Consider pairs of differential operators L{x, dj:),A{z, d^) in two variables x and 
for which there exists a function z) that is simultaneously a parametric family of 
eigenfunctions of L and A. 

LiIj{x,z) = f{z)ip{x,z) (1.1a) 
AiIj{x, z) = (j){x)'ijj{x, z), (l-lt>) 

where the eigenvalues f{z) and 0(x) are nonconstant functions of the variables z and 
X. Such pairs of operators were named bispectral and studied by Duistermaat and 
Griinbaum in [DG]. Choosing a normalization in which both L and A have unit lead- 
ing coefficients and constant next to leading coeffients (which is always possible, up 
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to reparametrization) , they were able to draw a number of conclusions about these 
operators. In particular, the coefficents of L and A are rational in the variables x and 
z, respectively, and the functions f{x) and 4>{^) ^ire polynomials. Furthermore, for the 
case of second order operators: 

L=^+.(x) (1.2) 

[DG] were able to determine all possible u{x); namely (up to translations in x or 
addition of a constant), 

c 

u{x) = X (Airy) or u{x) = — {Bessel) (1.3) 

or anything that can be obtained from the two cases 

uix) = 0, uix) = -^ (1.4) 

through the application of rational Darboux transformations. 

1.2 Bispectral Algebras of Rank 1. 

Wilson [Wl, W2] considered the case when L is embedded in a commutative 
algebra A of differential operators, sharing the same family of eigenfunctions 'ip{x,z), 
such that the orders of the elements of A are relatively prime (a rank 1 algebra) . Such 
algebras can be characterized by the associated spectral data, consisting of an algebraic 
curve, the spectral curve of A, denoted spec(.A) and, in general, a line bundle whose 
fibres are the joint eigenvectors. This data may be viewed as determining a point 
G Gr in the Hilbert space Grassmannian of Sato [SS] and Segal-Wilson [SW], such 
that the bispectral wave function ^(x, z) is the corresponding Baker- Akhiezer function 
tfjwi^, z) at the point t :— {ti,t2, ■ ■ ■) with (ti = x,ti = 0, i > 2). We recall that the 
Baker- Akhiezer function is the unique function of the form 

V^w^(t,z) = 7(t)(|l + ^^j, (1.10) 

where 7(t) := e^j=i ■'^ , taking values in the subspace W of the Hilbert space H := 
L'^{S^,C) of square integrable functions on the unit circle in the complex z-plane. 
Here t = {ti,t2, ■ ■ ■) is an infinite component vector and 7(t) may be viewed as an 
element of the infinite abelian group r_|_ consisting of elements of L^(S'^, C) admitting 
a holomorphic continuation to the interior of in the complex 2;-plane and taking 
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value 1 at the origin. The Baker-Akhiezer function may, in turn, be expressed in terms 
of the Sato tau function associated to W, by the formula [DJKM] 

Tw{t) 

where the components of [z] = {zi, Z2, ■ ■ ■) are 

Geometrically, tw is understood as the determinant, defined up to normalization, of 
the projection P_|_ : — > 7^+, where is the image under 7 G r_|_ of the subspace 
W C H and H+ C H is the subspace consisting of elements admitting a holomorphic 
extension to the enterior of 5"^. 

Wilson [Wl] gave the following equivalent characterizations of such VF's corre- 
sponding to bispectral algebras of rank 1. 

Theorem (Wilson [Wl]). The algebra A is bispectral and of rank 1 if and only if 
the following equivalent conditions hold: 

(i) spec{A) is a rational algebraic curve with only cusp singularities. 

(a) Within the proper normalization, the Baker-Akhiezer function 

ij{x, z) = iJw{x = tutj = 0,j> 2, z) = e"" (^^ + fl ^ j (l-^^) 

belongs to a plane W belonging to a subvariety Gr^"^ of the rational Segal-Wilson 
Grassmannian Gf^^, called the adelic Grassmannian (defined below). 

The adelic Grassmannian Gr^'^ is shown in [Wl] to consist of the PF-planes in the 
Hilbert space V. of the form 

W=-^Wl (1.14) 
q{z) 

where q{z) is a polynomial of suitable degree, rendering the projection P+ : — > H.^ 
to the subspace 7i+ (ZH a, Fredholm operator of index 0, and the subspace C 71+ 
is the intersection of the kernels of a finite number of linear forms {c^ G 7i!!j_, }, each of 
which has support at one point in the complex ^-plane. That is, each determines 
a condition of the form 



cfi9) = j:c^a9^''\Xi) = 0, (1.15) 



0=1 
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where g^°'^ denotes the a-th derivative of ^ G and {c^^} is a finite set of coefficients. 
The Aj's are just the roots of the polynomial q{z), with multiplicities coinciding with 
the number of conditions localized at that root. 

Wilson also showed that the property of bispectrality could be viewed as a con- 
sequence of the existence of a "bispectral involution" 

b : Gr^'* ^ Gr^"^ 

(1.16) 

b:W^W' 

such that 

'ipw'{x,z) = i/^wiz.x) (1.17) 

The larger, rational Grassmannian Gr''^* contains all rational and soliton solutions 
to the KP hierarchy. It was shown by Krichever [Kr], following earlier work of Airault, 
Mckean and Moser [AMM] on rational solutions to the KdV equation, that rational 
solutions of the KP equation which tend to zero as a; — > cxo can be expressed as 

where the location of the poles Xi{y,t) is determined by the fact that they satisfy 
the equations of the first two commuting flows of the Calogero-Moser system [M]. 
It follows that these are associated to bispectral tau functions rw- In [W2], Wilson 
showed that, by suitably completing the complexified Calogero-Moser phase space so as 
to allow for particle collisions, this extends to a complete parametrization of the adelic 
Grassmannian Gr*^*^, and hence the set of rank 1 bispectral algebras. More precisely, 
let Cn denote the completed n-particle rational Calogero-Moser phase space, defined 
as the space of pairs {X, Z) of complex n x n matrices whose commutator [Z, X] is a 
rank 1 perturbation of the identity matrix I, quotiented by the natural action of the 
general linear group 0[(n, C) on such pairs; 

g:{X,Z)^{gXg-\gZg-'), ge<Sl{n,C). (1.19) 

Define a map from the union of the spaces Cn to the adelic Grassmannian 

W : ^ Gr-d 

1.20) 

W:[{X,Z)]^W{X,Z), 
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where W{X, Z) e Gr^*^ is determined by the following formula for the corresponding 
Baker function at t = (x, 0, . . . ). 

^w{x,z) {x, z) = e"^ det (l - (X + xI)-\Z + ^I)-^) . (1.21) 

(Note that tpwixiz) and its derivatives at a: = span W.) Equivalently, the corre- 
sponding tau function may be expressed 

CO 

rwix,z) = det(X + xl + Y.M-Zy^)- (1-22) 

Wilson showed in [W2] that this identification gives an isomorphism of affine, non- 
singular, irreducible algebraic varieties. It is clear from (1.16), (1.21) that under this 
identification, the bispectral involution is equivalent to the symplectic involution 

h : Cn Cn 
b : {X,Z) ^ (Z^X*) 

on the completed Caloger-Moser phase space. (This was noted for the case of the 
original real Calogero-Moser phase by Kasman [K] , with b viewed as Moser's linearizing 
involution [M].) 

This therefore provides a complete determination of all wave functions ipw corre- 
sponding to rank 1 bispectral algebras. The following are two of the simplest examples 
of these. (In each case, the algebra A is identified with the space of polynomials which, 
under multiplication, leave W{X, Z) invariant. Only the simplest representative ele- 
ments L and A, are given.) 

Example (i) 



(1.23) 



n=l, X={a), Z={0) 

1 



^w{x,z){x,z) 
L 
A 



{x + a)z 
2 



dx"^ {x + a)"^ 
^ ^ d 
dz^ 



+ 2a- 7T, (j)(x)—x + 2ax 

dz z^ 



(1.24a) 
(1.24b) 

(1.24c) 

(1.24d) 
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Example (ii) 



d^ 6 d 12 ^ n ... ,^ 



2. Isomonodromic Deformations. 

2.1 Rational Covariant Derivatives Operators. 

Consider a parametric family of matrix first order differential operators 

where B{X,r) is an r x r matrix valued polynomial in A that may depend also on a de- 
formation parameter r, {Nia}i<i<n,i<a<mi is a set of r x r matrices, also dependent on 
T, and {ai{r) G C}i<i<n are the location of the finite poles, also possibly r dependent. 
This may be viewed as a r-parametric family of rational covariant derivative operators 
defined on a rank r vector bundle over the Riemann sphere, punctured at the poles 
(including possibly oo). The generalized monodromy data associated to the operator 
T>x consists of the local monodromy matrices about each of the poles, together with 
the Stokes matrices and connection matrices [JMU] . Assuming local differentiability in 
the deformation parameter r, it is a basic fact that this monodromy data is invariant 
under variations in r provided there exists a second differential operator 

V,:= — -E{X,t), (2.2) 

where -^(A, r) is also a r-parametric family of r x r matrices, rational in A, such the 
commutativity condition 

[Vr,Vx] = (2.3) 

is satisfied. This means that the data {B{X, r), af(r), Niair)} determining T>\{t) satisfy 
a set of first order ODE's in the deformation parameter r. It also implies that, locally 
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at least, an invertible matrix valued function ^{X,t) exists, uniquely determined by 
any given set of initial values, simultaneously satisfying 

Pa* = (2.4a) 
P^* = 0. (2.4b) 



The simplest case, studied since the beginning of this century, is when the operator 
T>x is Fuchsian; i.e., it has only regular singular points. This means that all the finite 
poles are simple (mj = 1, Vz) and B{X, r) = 0, so T>\ is of the form 

OA ^-^ A — Ctj 

1=1 

For this case, it has long been known that the only way the residue matrices {-^i}i=i,...n 
can vary is through their dependence on the pole locations {Q:i}i=i,...n: and they are 
constrained to satisfy an overdetermined system of PDE's. Choosing the a^'s as the 
deformation parameters, the associated infinitesimal isomonodromic deformation oper- 
ators are 

f..:=^ + v^. (2.6) 

The corresponding commutativity conditions 

[PA,I^aJ=0, [D«,,D„J=0, l<,i,j<n (2.7) 

are mutually compatible (i.e. integrable in the Frobenius sense), and equivalent to the 
following system of equations, known as the Schlesinger equations 



dN, _ [iV„iV,] 



i+) (2.8a) 



^ = (2.8b) 

(In particular, the simplest Fuchsian case involving nontrivial isomonodromic deforma- 
tions is when n = 3 and r = 2, for which these reduce to the equations of the sixth 
Painleve transcendent Pyi-) 

The corresponding determination of all equations of the type (2.3) generating 
isomonodromic deformations of operators of the general form (2.1) was made in [JMU]. 
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For the simplest nonfuchsian case, where the finite poles remain first order, but an 
irregular singularity of Poincare index 1 is added at infinity by allowing S to be a 
nonzero diagonal matrix B{t) = diag (/?i(t), . . ./^^(t)) with distinct eigenvalues that 
are independent of A, we have 

Vx{r) = -^-Af{X,T) (2.9) 

where 

A^(A,r):=S(r) + X:^. (2-10) 

In this CctSGj ctS shown in [JMMS], besides the original ctj's, the possible deformations 
parameters include the eigenvalues {Pa}i<a<r of the matrix B. The corresponding 
infinitesimal isomonodromic deformation operators are 

2^a,:=^ + i-^, ^ = l,...,r^ (2.11a) 
oai A — ai 

OPa Pa Pb (2.11b) 

where Ea is the elementary matrix with 1 in the aa position and O's elsewhere, and 

n 

N^:=J2Ni. (2.12) 

i=l 

Again the mutual commutativity of the operators {X'a{)^^/3a)^^A}i=i,...n,a=i,...r) 

[V^„ V^.] = [Vo^,M = [V^^, Vp^] = [PA,I^aJ = pA,I^/3j = 0, (2.13) 

gives a Probenius intcgrable system for the residue matrices Ni, locally determining 
their dependence on the parameters {ai,(5a\- 

2.2 R-Matrix theory and Duality. 

The above isomonodromic deformation equations may be viewed as nonauto- 
nomous Hamiltonian equations ([JM], [H]) on the space (0l(r, C))"^ of n-tuples of r x r 
matrices {A^i, . . . N^}, with respect to the Lie Poisson bracket structure 

{{Ni)ab, {Nj),d} = Sij [SabNcb - SbcNda] ■ (2.14) 
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Equivalently, the rational matrix valued functions M{X), for different values of the 
complex parameters A and fj,, satisfy 

{J\fiX)Wil^)} = [r(A - XiX) ® I + 1 Xifi)] , (2.15) 

where r(A — /x) is the rational i?-matrix, viewed as an element of End(C"'0C"'), defined 

by 

r(A) := (2.16) 

where 

Pi2(u(8)v) = v®u (2.17) 

is the transpostion endomorphism. This defines a Poisson bracket structure on the space 
g[^^^(r, C) of rational r x r matrix valued functions of A. The a^'s and /5a's may be 
viewed as multi-time parameters associated to the n + r nonautonomous Hamiltonians 
{Hi, Ka} i=i, - ri defined by 

a=l, .. .r 

H,:= — (( tr (Af\\)) d\ = tr{BNi) + V (2.18a) 

K^:=J-I dXxi dz\tr{{B-zl)-^M{X)f -ztT{{B-zl)-^M{X)) 

= E -,{N,U + E ^""-^^'^y^^ (2.18b) 



Defining a 1-form 9 on the parameter space by 

n r 

d:=J2 Hidai + E ^adf^a, (2.19) 

the following results, which may be verified directly [JMU], are a consequence of the 
general i2-matrix theory, adapted to the nonautonomous setting [H]. 

Theorem. The isomonodromic deformation equations (2.13) are Hamilton's equations 
for the set of Hamiltonians {Hi, Kg} contained in the 1-form 9 

■■= E + i2dPa^ = {Ni, 9}. (2.20) 
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Furthermore, the Hamiltonians {Hi, Ka} all Poisson commute amongst themselves, and 
hence 6 is a closed form on the parameter space 

de = 0, (2.21) 

and there exists, at least locally, a function r (the tau function) such that 

d = d{lnT). (2.22) 



The Hamiltonian equations determining these isomonodromic deformations may 
in fact be hfted to an associated symplectic vector space, and then projected to another 
space of rational differential operators, giving rise to a second, dual representation of 
these equations as isomonodromic deformations [H]. It is easy to see that all j\/'(A)'s of 
the form (2.10) may be expressed as 

J\f{\) ^ B + G^{A- Xiy^F, (2.23) 

where ^ is a diagonal N x N matrix, 

A:=diag (ai,...,an), (2.24) 

with N = Yl^=i — rk(-/Vi), whose eigenvalues {ai\ are at the poles of and 

have multiplicity ki, and (F, G) are a pair of rectangular N x r matrices. If the space 
j^Nxr ._ ^^p^ Q-^y such pairs is given the canonical symplectic structure 

uj = ti{dF^ A dG), (2.25) 

the map : M.^^"^ 0^rat(^; 'the space 0trat(^; defined, for fixed A and B, 

by 

J^:{F,G)^AfiX), (2.26) 

is a Poisson map. In fact, it may be viewed as defining a Poisson quotient of A4^^'^ by 
the canonical Hamiltonian action of the stability subgroup &a C &1{N, C) of A, given 

by 

{g,{F,G))^{gF,ig^)-'G). (2.27) 



(The ki X r blocks {Fi,Gi) within the matrices F and G corresponding to the eigen- 
values ctj determine the local monodromy of the operator Vx about these points.) 
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Denoting the pull-back of the 1-form 6 under by ^, the isomononodromic defor- 
mation equations (2.13) may be viewed as the projection to Qlratir^ C) of the following 
nonautonomous Hamiltonian equations on the auxiliary space M.^^'^: 

dF = {F, 6}, dG= {G, 6}, (2.28) 

where 

d = T.d^^a^ + J2dPa^. (2.29) 

j=l a=l 

To obtain the dual set of isomonodromic deformation equations, we define a new 
Hamiltonian quotient of M.^^^^ this time by the stabilizer 25^ C S5l(r, C) of the 
element B e S5l(r), by the dual Poisson map : (F, G) — > M.{z), where 

''' nf 

M{z) :=A + F{B - zlrY^G'^ = ^ + J] (2-30) 

Defining the associated set of differential operators 

V,:=^^-M{z) (2.31a) 

ooLi ^ cti - aj 

Hamilton's equations (2.28) are projectible under the map and imply the commu- 
tativity conditions 

V^,] = [Da^^fiJ = [Va^, VpJ = [D^.V^,] = [D^VpJ = 0, (2.32) 

which, in turn, imply the invariance of the monodromy data of the operator "D^ (resp. 
Vx) under changes inthe parameter x (resp. z). These equations may also be inter- 
preted as Hamiltonian equations on the space of rational NxN matrix valued functions 
of A4 (z) , generated by the set of Poisson commuting, nonautonomous Hamiltonians ob- 
tained by projection: 

Ka:=^.i tr [M\z)) dz = tr(Ai?„) + ^-^f^ (2.33a) 

dzz (f dX [tr {{A - XlNr^M{z))^ - Atr {{A - XIn)~^M{z)) 

= Tf3b{M,h + T (2.33b) 
^-^ ^-^ a4 — a ■ 



6=1 3=1 



■3 
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The duality map 



(F, G, A, B, A, z) ^ (G^, F^, B, A, z, A) 



(2.34) 



may be viewed, after projection, as transforming the system of isomonodromic defor- 
mations equations (2.13) into the dual system (2.32). 

Tthis duality map is reminiscent of the bispectral involution (1.23) expressed in 
terms of the matrices (X, Z) determining the Calogero-Moser system. A full explana- 
tion of the relationship between these two maps is not yet developed, but one certainly 
exists. What will be provided in the next section will be a reformulation of the two 
examples of bispectral pairs given in section 1 in terms of an equivalent, dual pair of 
isomonodromic deformation equations - in the very particular case where the mon- 
odromy happens to be trivial. These examples in fact illustrate a correspondence that 
can be made between the bispectral pairs of the type studied by Wilson; i.e., thoe 
belonging to rank 1 bispectral algebras, and dual pairs of parametric families of ma- 
trix differential operators of first order, depending rationally on both parameters, and 
having trivial monodromy. 

3. Relation between Bispectrality and Dual Isomonodromy. 

The general structure of bispectral equations and infinitesimal isomonodromic 
deformation equations, both involving a wave function depending on two variables 
that simultaneously satisfies a linear differential equation in each of them, suggests 
that the two may be related. The similarity between the duality map (2.34) and the 
bisectral involution (1.23) further suggests a relation between bispectrality and duality. 
The deeper relation between bispectral pairs and dual isomonodromic deformations 
involves the theory of dressing transformations and the tau function, and will be not 
developed here. However, to demonstrate that such a relation does exist, we consider 
the two examples of bispectral rank 1 operators introduced in section 1.2, and show 
how an equivalent infinitesimal isomonodromic deformation system can be derived for 
each. 

3.1 Examples. 

(i) For the case of the Baker type bispectral wave function 



there exists a second bispectral wave function for the same pair of operators L and A 
given in (1.24c,d), namely 




(3.1) 



■01 := '4'{—x — 2a, z) = e 



-2az 



ip{x,-z) 



(3.2) 
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Forming the Wronskian matrix 



we may associate the following pair of rational 2x2 covariant derivative operators 

5/0 



d / -a \ 

^'^=S-((a + x).+ ^ r-a). (3.4b) 

which simultaneously annihilate ^, 

V,^ = 0, V,^ = 0, (3.5) 

This is equivalent to the bispectral conditions (l.la,b) for the operators L, A defined 
in (1.24c,d) and implies the commutativity condition 

[D^, V,] = 0. (3.6) 

It follows that the monodromy of either of these operators, say T>g, is invariant under 
the deformations determined by varying the other parameter, x. Note however the fact 
that there is not just one single-valued bispectral wave function, but a basis of single 
valued bispectral wave functions, which implies in this case that the operators Vx and 
actually have trivial monodromy for all parameter values. 

To obtain the dual system of (trivial) isomonodromic deformation equations, we 
just define new operators similarly, using the Wronskian matrix formed by differenti- 
ating with respect to the variable z: 

This is simultaneously annihilated by the operators 

d 



\ a+x {a.+x)z a+x 
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implying tlie commutativity condition 



[Da^, V,] = 0, 



(3.9) 



from which again follows that the monodromy of either of these operators is invariant 
under the deformations determined by varying the other parameter. 

Note however that the singularity structure of the operators , , and 
is such that the Hamiltonian i?-matrix formulation discussed in section 2.2 does not 
suffice. In order to include these, the i?-matrix approach must be extended (cf. [HTW] , 
[HR], [HWl], [HW2]) to allow for singularities of higher order at oo. 

(ii) For the case of the bispectral wave function 



because of the invariance of the eigenvector equations for the operators L and A in 
(1.25c,d) under the cyclic group generated by 




(3.10) 



Z I— > U)Z, u> 



:= e 3 



(3.11) 



we may define a basis {i/jq, V'l, V'2} of bispectral wave functions by 



'4)j{x,z) :^ip{x,u;^z), 0,1,2. 



(3.12) 



The Wronskian matrix 




^2 



(3.13) 



is then simultaneously annihilated by the operators 




(3.14b) 



(3.14a) 



This is equivalent to the eigenvector equations for the operators L and A of eqs. 
(1.25c,d), and implies the commutativity of the operators Vx,T>z. Hence the mon- 
odromy of each of these operators (which again is trivial because of the single valuedness 
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of the matrix "^), remains invariant under changes in the parameter values. Because 
of the symmetry of the bispectral wave functions under the interchange x z, the 
corresponding "dual" isomonodromic system is obtained by simply interchanging the 
two variables in the definitions of the operators Vx, T>g. As in the previous example, 
the proper treatment of these systems within the i2-matrix framework requires the 
inclusion of higher order singularities at oo. 

Proceeding similarly, we may associate to every rank 1 bispectral pair two mutu- 
ally dual 1-parameter families of differential operators having trivial monodromy for 
all parameter values, since in each case we may form a basis of single- valued bispectral 
wave functions. This procedure may also be applied to the case of higher rank bispec- 
tral operators, except that the monodromy of the resulting dual families of operators 
is no longer necessarily trivial, and wc obtain instead 1-parameter families of operators 
with constant monodromy under variation of the deformation parameters. These more 
general cases, as well as their relations to the evaluation of certain limits of Fredholm 
determinants, will be dealt with elsewhere. 
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